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INTRODUCTION 


Deardorff's  (1972)  similarity  theory  for  convective  boundary 
layers  (CBL)s  suggests  that  CBL  parameters  are  controlled  by 
buoyancy,  ft  =  g/0,  (where  g  is  gravitational  acceleration  and  6 
is  potential  temperature),  surface  temperature  flux,  w'e'0, 
height  z  above  the  surface,  and  zif  the  height  of  the  inversion 
base  which  marks  the  top  of  the  mixing-layer  (ML).  The  velocity, 
temperature,  and  humidity  scales:  w,,  u,ML,  0,N  ,  q„N1L,  formed  from 
these  parameters  are  used  in  universal  relationships  for  the 
gradients  and  fluxes  of  mean  and  turbulent  quantities  in  the  CBL. 
The  first  of  these,  the  convective  turbulence  scaling  velocity 


w,  =  (ftz,  w'e '0) 13  ,  (i) 

determines  much  of  CBL  physics,  including  diffusion.  Along  with 
w'0'o,  w,  forms  the  basis  for  the  remaining  parameters  used  to 
scale  the  heat,  moisture,  and  momentum,  as  well  as  passive  scalar 
fluxes.  However,  w,  neglects  mechanical  turbulence  due  to  shear 
and  is  only  strictly  valid  for  free  convection  induced  by  surface 
heating.  Yet,  w»  is  used  often  in  mixed  forced/free  convection 
cases  for  lack  of  alternatives.  Here  we  seek  w„s,  an  amendment 
to  w. ,  which  includes  mechanical  turbulence  contributions  under 
mixed  forced/free  convection  conditions.  We  assume  mechanical 
generation  comes  from  three  sources:  the  barotropic,  horizontally 
homogeneous  surface  layer  (SL) ,  barotropic  boundary  layer 
entrainment  zone  (EZ) ,  and  a  baroclinity  (BC)  component.  Beyond 
the  z  and  z,  length  scales,  we  propose  that  mixed  forced/free 
convection  also  depends  on  the  SL  roughness  length,  z0,  SL 
Obukhov  length,  L  =  -u,30/ (gkw 1 0 '0)  (where  u,  is  friction  velocity 
and  k  =  0.4  is  the  von  Karmen  constant),  EZ  thickness,  Ah,  and 
zh,  a  baroclinic  length  scale  which  we  define  below.  We  model 
these  mechanical  sources,  using  relevant  ratios  of  the  length 
scales,  thereby  fulfilling  similarity  requirements. 


APPROACH 

We  outline  the  basic  approach  before  discussing  details.  For 
non-baroclinic,  horizontally  homogeneous,  steady  state  cases,  the 
vertical  profiles  of  the  buoyancy  and  momentum  fluxes  must  be 
linear,  if  the  shapes  of  the  vertical  profiles  of  potential 
temperature  and  wind  speed  are  to  remain  static.  That  is, 


ftw '  0 '  /  -  ftw '  0 '0  ( 1  -  (1  +  a,)z/Zj)  ,  (2a) 


-u.'w*,  =  -u,'w'0(l  -  (1  +  a;)z/z,)  ,  (2b) 
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where  a[  =  -w '  0  '  zi/w '  0 '0  and  a2  =  u/w'^/Uj'w'o  (see  fig.  1).  The 
subscripts  "zi"  and  "0"  refer  to  conditions  at  the  inversion  base 
and  surface.  -u'w'0  =  u,:,  w'0'zi  -  WeA0,  and  -Uj'w'^  =  iU,We, 
where  the  CBL  entrainment  rate,  We  =  zit  and  a0  is  the  potential 
temperature  jump  across  the  EZ .  The  subscript  ",t"  denotes  the 
partial  derivative  with  respect  to  time.  Such  notation  will 
refer  to  spatial  and  temporal  derivatives  in  general.  The  Betts 
(1974)  correction,  w'0'zi  =  (a0  -  ruAh)We,  hardly  affects  eqn.  2a 
since  Kamada  (1988a)  showed  that  the  much  smaller  interfacial 
thickness,  lm,  must  replace  the  EZ  depth,  Ah.  For  lack  of  better 
information,  we  assume  that  above  z{  the  barotropic  flux  profiles 
are  linear  with  height  up  to  H  =  zx  +  Ah,  the  top  of  the  boundary 
layer.  I . e. , 

Z  “  z 

Bw  1  0  '  z  =  Bw  '  0  '  Z1  ( 1  -  - -  )  (3a) 

Ah 


Z  “  z 

-U.'w  'z  =  -U,'W'ZI(1  -  -  )  .  (3b) 

Ah 


Using  eqns.  (1),  (2a),  and  (3a),  the  vertically  integrated 

buoyancy  flux  or  total  CBL  buoyancy  production  rate  of  turbulence 
kinetic  energy  (TKE)  can  be  written  to  show  that  the  purely  free 
convective  turbulent  scaling  velocity  is  proportional  to  the  cube 
root  of  the  boundary  layer  integrated  TKE  production  rate, 

h  _  1  -  a,  ( 1  +  a) 

f  6w'0'z  dz  =  n  w,3  ,  where  n  =  -  , 

z0  2 

(4) 

and  a  =  Ah/z,.  0  =  0.4  for  mature  CBLs  because  under  such 
conditions,  a,  -  0.2  and  a  ~  0.1.  Since  this  may  not  be  true  for 
other  conditions,  we  discuss  D  in  detail  below.  However,  here  we 
propose  that  this  sort  of  proportionality  extends  to  shear;  that 
mechanical  contributions  to  a  shear  inclusive  turbulent  rraling 
velocity,  w.s,  may  be  written  in  terms  of  the  boundary  1  jyer 
integrated  mechanical  generation  rate  of  TKE.  I.e., 


w 


*s 


1  1/3 

H  _  H _ 

j  Bw ' 0  ' ,  dz  +  j-u1iw'ui  z5z 
z0  zO 

n 


(5) 
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fig.  1  Schematic  profiles  of  a)  buoyancy,  b)  stress,  c)  U  velocity, 
and  d)  TKFi  shear  generation  rate  for  barotropic  case 


-u'w'  L’  "  -u'w’dU  dz 


We  also  propose  that  the  CBL  shear  generation  rate  can  be 
represented  by  discrete  barotropic  surface  layer  (SL) , 
entrainment  zone  (EZ) ,  and  baroclinic  (BC)  components. 


SL  EZ  BC 

H _  -L _  H  _  H  _ 

f-Uj'w'u^Sz  =  f  -u, '  w '  U|  Zdz  +  |-ui’w'uiz3z  +  j-u, 'w'ulb  Zdz 

zO  zO  zi  0 

(6) 

Here  uib  z  is  just  the  baroclinic  portion  of  the  vertical  shear. 

BAROTROPIC  SURFACE  LAYER  SHEAR 

We  first  consider  the  barotropic  SL  and  assume  that  most  of  the 
shear  occurs  below  the  height,  -L.  From  SL  similarity  we  have 

u  z  =  (u./kz )0m  ,  with  0m  =  (1  -  15z/L)13  ,  (7a,  b) 

(Carl  et  al.  1973).  Some  workers  prefer  0m  formulas  involving 
the  -1/4  power.  However,  the  choice  is  among  empirical 
expressions  and  integrals  of  -1/4  power  expressions  result  in 
awkward  complex  numbers.  The  accepted  value  in  the  free 
convection  limit  is  -1/3,  so  we  choose  this  for  theoretical 
appeal  as  well  as  mathematical  convenience. 

Since  <pm  is  only  valid  for  z  >>  z0,  we  set  the  upper  integration 
limit  at  -L,  the  lower  limit  at  a  somewhat  arbitrary  10  z0,  and 
treat  the  layer  below  10  z0  separately.  Above  10  z0,  from  eqns. 
(2b)  and  (7a),  we  have 


-L _ 

\  -Ul’w'uizaz 


10  zO 


u,3  -l  0ma  z 

—  J  — 

k  lOzO  z 


u«(u»:  -  AU,We)  -l 

-  1  0maz 

kZj  lOzO 


U,3  -L  0maz  -L 

- \  -  -  (l-a;)  (L/z,;)  w,3  l  <pmdz 

k  10z0  Z  10z0  .  (8) 


From  the  L  definition,  u,3/k  =  (-L/Zi)w,3.  With  the  axes  aligned 
with  the  mean  wind,  we  have  -u'jW'uz  =  -u'w'uz,  since  (vz,  wz)  << 
uz  in  the  mixed  layer.  If  shear  supplies  sufficient  TKE,  10  z0 
<<<  -L.  So  in  eqn.  (8)  term  (2)  =  0.54  (1  -  a;)  (L/z,)2w,3. 

Letting  p  =  -I5z/L,  the  first  integral  reduces  to 
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V3  (2  (p+1)  1/3  +  1) 


Ln [  (p+l)2'3  +  (p+1)  1/3+  1] 


=  V3  ATAN 


3 


-  + 

2 


Ln [  (p+1)  1/3 


1] 


-L 

lOzO 


The  upper  limit  evaluates  to  =  1.51.  The  first  two  terms  in  the 
lower  limit  =  1.26,  if  L  >>  150z0.  So  in  eqn.  (8|  the  first 
integral  =  0.25  -  Ln  R,  where  R  =  (1  -  150z0/L)  1  3  -  1  <<  1.  A 
5cm  z0  and  150m  L  would  give  150  z0  =  -0.05L.  Above  10  z0,  this 
leads  to  =  4.36u,3/k  from  term  (1)  in  eqn.  (8).  Oddly,  R  varies 
inversely  with  roughness  length.  But  from  u,  =  ku/(ln(z/z0  -  <^m) 
and  the  definition  of  L,  we  see  that 


nz0/L  =  (ngw'0' /k2u3©)  z0(ln(z/z0)  -  <^m)3 


(9) 


So  for  given  levels  of  u,  w'Q',  0,  and  z,  the  ratio  nz0/L  and 
thus  term  (1)  may  vary  much  less  than  z0  itself. 

Below  10  z0  we  neglect  non-neutral  contributions  to  u.  and  assume 
a  constant  value,  -u/w'u,  =  u,3/10kz0.  Then 

10z0  _ 

j  -u,'w'uzdz  =  u,3/k  .  (10) 

o 

Together,  eqns .  (8)  and  (10)  give  an  estimate  of  the  vertically 

integrated  SL  shear  production  rate, 


j  -ut'w'uzdz  =  (1.25  -  Ln  R)u,3/k  , 

o 

=  (-w,3L/z,)  (1.25  -  Ln  R)  .  (lla,b) 


Note  that  we  have  dropped  term  (2)  from  eqn.  (8),  since  1)  the 
coefficient  is  an  order  of  magnitude  smaller  than  for  term  (1), 
2)  term  (2)  is  second  order  in  z;,  3)  usually,  -L  <  0.4  z,,  and 
4)  a:  is  of  order  10'1  -  1;  so  some  further  cancellation  occurs. 

From  eqns.  (2)  and  (3),  if  fl  =  0.4,  we  have  the  following  ratio 
of  enhanced  to  unenhanced  TKE  contributions, 
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1 1/3 

(1.25  -  Ln  R)L/Zj 

-  w,  =  (1  +  (3.1  -  2. 5Ln  R)  (-L/z,)  )  !  -  w.  . 

n 

J  (12) 

For  our  previous  example,  w,s  exceeds  w,  by  a  factor  of  =  (1  + 

3 3u.’/w.3) 1/3 .  In  this  case,  since  u*/w.  =  0.54,  SL  shear 
generation  alone  may  enhance  w,s  by  more  than  80%. 

To  be  useful,  w»s  should  maintain  the  same  proportion  to  the  mean 
turbulent  velocity,  auj/3,  in  both  the  convective  and  neutral 
limits.  We  test  this  by  noting  for  the  barotropic  case  that  the 
cube  root  of  the  vertically  integrated  production  rate  is  usually 
—  2.4  u,  for  SL  shear.  The  EZ  and  BC  terms  are  not  relevant  for 
the  pure  neutral  case.  From  Grant  (1986)  and  Mason  and  Thompson 
(1987),  we  find  that  aui/ 3  =  1.3  u,  in  the  layer  averaged  neutral 
boundary  layer.  This  implies  that  aui/3  =  0.53  w»s.  In  the  pure 
convective  limit,  the  CBL  averaged  value  for  aui/3  is  ~  0.55  w, 
(Caughey  and  Palmer,  1979  and  Deardorff  et  al.  1980).  So,  the 
ratio  of  w«s  to  the  layer  averaged  mean  turbulent  velocity  is  the 
same  in  the  convective  and  neutral  limits,  as  required. 


ENTRAINMENT  ZONE  SHEAR 

The  EZ  also  generates  mechanical  turbulence.  To  quantify  it  we 
assume,  moreover,  that  shear  is  constant  across  Ah.  Thus, 

h  _ 

j  -Uj'w'uitZdz  =  WeAUf/2  .  (13) 

zi 

Equation  (13)  suggests  the  need  for  a  simple  shear  inclusive 
model  to  estimate  the  entrainment  rate,  We.  For  this  purpose  we 
parameterize  the  TKE  equation, 


storage  buoyancy  transport  shear  dissipation 

e(  =  ft  w'8'  +  (u^e  +  p  'e/p)  z  -  Uj'Uj'Ujj 

at  the  inversion  base  as, 


storage  buoyancy  transport  shear  dissipation 

Cs!w,;We/Ah  =  -ftA0We  +  Ccw,3/z,  +  CmAU,:We/Ah  -  Cdw,?/z, 

(15) 
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The  scaling  arguments  are  that  1)  again  w'9'zi  =  -WeA0  and  2) 
large  mixed  layer  eddies  drive  the  TKE  transport  rate  into  the 
EZ.  Thus  TKE  transport  should  vary  with  the  TKE  (a  w«2)  times 
the  large  eddy  turnover  time  scale,  zj w,.  3)  The  shear 
production  rate  is  modeled  as  the  momentum  flux,  AU,We,  tines  the 
local  EZ  shear,  AU,/Ah.  4)  To  maintain  a  steady  state  inertial 
subrange,  the  dissipation  plus  buoyancy  destruction  rates  must 
balance  the  large  eddy  TKE  injection  rate  due  to  transport  and 
shear  production.  So  in  the  neutral  limit,  dissipation  should 
scale  like  transport  with  the  change  rate  of  large  eddy  TKE, 
w,3/Zj.  Since  dissipation  occurs  at  small  scales,  unaffected  by 
buoyancy  anisotropy,  we  retain  the  neutral  form.  For  simplicity 
we  neglect  the  small  TKE  storage  rate  term.  Cm  estimates  vary 
somewhat,  (Pollard  et  al.  1981;  Kato  and  Phillips,  1967,  Tennekes 
and  Driedonks,  1981),  but  we  see  little  reason  for  departures  far 
from  unity,  so  we  assume  Cm  ~  1.0.  We  discuss  Cc  and  Cd  below. 

Tennekes  and  Driedonks  (1981)  presented  a  similar  model,  except 
that  here  we  scale  the  shear  production  locally  with  Ah  instead 
of  z,.  Solving  for  the  entrainment  rate  gives, 


We 


(Cc  -  Cd) 

- - -  w,3/z, 

Ba©  -  Auy/Ah 


(16) 


Therefore,  from  eqn.  (13)  we  have 


H  _ 

i  -u/w'u.^z 

zi 


(Cc  -  Cj)  Ah 

-  -  w,3 

2(Rb  -  1)  z, 


(17) 


where  the  bulk  Richardson  number, 

Ba  0a  h 

Rb  =  - —  .  (18) 

au," 

The  EZ  is  transitionally  turbulent  by  nature.  Thus,  we  expect  Rh 
to  fluctuate  about  some  critical  threshold  turbulence  value,  Rch. 
Rcb  is  much  larger  than  the  1/4  value  accepted  for  the  stable  SL 
because  EZ  TKE  is  mainly  maintained  by  large  eddy  transport  from 
below.  So  shear  generation  is  only  part  of  the  TKE  injection 
rate,  and  the  dominant  eddies  are  also  much  larger  than  in  the 
transitionally  turbulent  SL.  Thus,  the  EZ's  large  negative 
buoyancy  will  destroy  much  of  the  TKE  during  the  eddy  cascade 
before  it  can  be  dissipated  at  small  scales.  For  these  reasons 
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Rcb  may  exceed  unity  in  the  EZ . 

Indeed,  our  simulations  to  steady  state  using  such  an  entrainment 
model  suggest  a  typical  value  of  Rcb  =  1.3.  From  the  results  of 
Turner  (1968),  Deardorff  et  al.  (1980),  and  many  others,  we 
suggest  that  Cc  -  Cd  =  0.25  in  the  convective  limit.  So  in  eqn. 
(17)  the  EZ  shear  contribution  =  0.5  aw,3,  where  a  s  Ah/z,. 

To  assume  that  a  —  0.1  -  0.2  may  often  suffice.  Then,  a!  =  0.2, 
n  =  0.4,  and  the  barotropic  EZ  contribution  to  w.s  is  only  a  few 
percent.  But  a  can  exceed  unity  during  rapid  entrainment  into  a 
near-neutral  layer  remaining  from  the  previous  day's  CBL  (Nelson 
et  al.  1989).  During  such  times,  entrainment  models  suggest  that 
a,  may  exceed  40%  (Kamada,  1988a).  Then  n  —  0.1,  and  EZ 
enhancement  of  w,s  may  exceed  0.4.  With  SL  shear  this  may  cause 
w,s  to  double.  In  this  case  a  better  diagnostic  or  prognostic 
estimate  for  a  is  needed.  For  a  diagnostic  estimate,  we  begin  by 
assuming  Rb  =  Rcb  and  solve  eqn.  (17)  for  Ah.  This  leaves  unknown 
au(  and  a©,  the  inversion  temperature  jump.  a©  is  hard  to 
specify  more  precisely  than  -  +/-  30%  because  Z;  and  H  are  not 
that  highly  resolved.  Yet,  we  know  that  a©  grows  with  the 
boundary  layer  growth  rate,  Ht,  into  the  layer  aloft  (with  lapse 
rate  ru)  ,  and  shrinks  with  the  mixed  layer  warming  rate,  ©m>l. 

So,  using  the  simplest  steady  state  assumption, 


a©  ,  = 


ruH,  - 


© 


m.t 


=  o 


(19) 


The  mixed  layer  warming  rate  is  given  by  the  surface  heating 
rate,  plus  the  entrainment  rate  of  warm  air  from  aloft,  i.e., 


©m.t  =  ( W  '  ©  *  o  -  W'Q'J/Z,  .  (20) 


So,  if  we  approximate  H(  as  We ( 1  +  a) ,  we  have 


a©  =  w'©'0/We  -  ruz,(l  +  a) 


(21) 


Using  eqns.  (16)  and  (18)  and  combining  constants,  this  results 
in  two  solutions  for  a.  We  show  only  the  physical  root, 


a 


2  .  3 au," 

ru  fl  z,: 

L 


+ 


4 


J 


1 


2 


(22) 
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Though  not  explicit,  eqn.  (22)  implicitly  involves  a©  because  it 
presumes  through  eqns.  (18  -  22)  that  a©  depends  on  Tu,  We, 
w'0'o,  and  Z;.  This  assumes  that  the  capping  inversion  requires 
finite  potential  temperature  and  windspeed  jumps,  whenever  there 
is  horizontal  wind. 

In  eqn.  (22)  we  must  still  estimate  the  windspeed  jumps  across 
the  EZ.  Wyngaard  (1988)  proposed  the  following  first  order 
estimate.  For  steady  flow  over  a  horizontally  homogeneous 
surface,  the  stress  gradient  is  given  by 


-u'w'  z  =  f(v  -  vg)  ,  and  -v'w'z  =  f  (ug  -  u) 


(23a, b) 


where  f  refers  to  Coriolis  forcing,  and  the  subscript,  g,  refers 
to  geostropic  winds  (Panofsky  and  Dutton,  1984).  For  barotropic 
flow  with  the  axes  aligned  with  the  mean  mixed-layer  wind, 


-u'w'  z  =  fvg  =  (We(ug  -  u)  +  u,:)/z1  ,  and 

-v'w'  z  =  f(ug  -  u)  =  We(vg  -  v)  /z,  .  (24a, b) 

The  momentum  balance  becomes 

We(ug  -  u)  +  Ut"  =  -fz ,vg  ,  and 


-We  Vg  =  f  z j  ( ug  -  u) 


(25a, b) 


For  barotropic  cases,  if  au  =  (ug  -  u)  and  av  =  (vg  -  v)  across 
the  EZ,  then  the  solutions  to  eqns.  (25a, b)  are  to  first  order, 


au  =  We(u,/fz,):  ,  and 


AV  =  U.'/fZ, 


(26a, b) 


Here,  almost  all  of  the  EZ  shear  results  from  turning  the  wind 
into  the  v  direction.  So  for  a  mature  mid-latitude  BL  with  u,  = 
0.4ms1,  f  =  10 V,  z,  =  10-m,  0  =  300  °C,  and  Tu  =  3  x  10’3  °C  m1, 
we  have  a  =  0.057.  However,  during  growth  into  the  near  neutral 
remnant  mentioned  earlier,  conditions  may  be  more  like  Tu  =  1  x 
1 0'3  °C  m'1,  and  z,  =  500m.  This  would  result  in  a  =  1.3.  Values 
spanning  such  a  range  were  reported  in  Nelson  et  al.  (1989). 

They  also  observed  and  prognostically  modeled  the  typical  early 
morning  to  late  afternoon  hysteresis  found  in  plots  of  a  versus 
We/w,.  They  noted  that  such  hysteresis  is  absent  from  previous 
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diagnostic  expressions  for  the  EZ  depth  which  have  the  form  a  <x 
Ri,n,  where  Ri,  e  fiz^e/w,2  and  1/4  <  n  <  1.  No  hysteresis  occurs 
in  such  expressions,  if  we  assume  (see  eqn.  2a)  that  a,,  the 
inversion  to  surface  heat  flux  ratio,  is  constant  because  then 
We/w,  oc  l/Ri,.  Unlike  such  expressions,  eqn.  (22)  shows 
qualitatively  appropriate  hysteresis.  In  fact,  by  combining 
eqns.  (16),  (21),  (22),  and  (26),  we  see  that 

a  u,4w, 

-  =  - ; -  .  (27) 

We/w,  w'6,06f2zl3 


Far  from  being  constant,  aw, /We  will  typically  increase  initially 
as  surface  heating  and  u,  grow  in  the  young  day-time  boundary 
layer,  but  will  diminish  later  as  z-  becomes  large. 

If  we  require  more  accuracy,  we  can  couple  eqns.  (16,  18,  19,  20, 
and  26)  with 


j  2u.:av  We  a  v2a  0  ,  j 

H  t  =  We  +  (R,b/J5)  | - ; - —  |  (28) 

j  f  Z  j”A  0  A02  j 

L  J 

to  provide  a  complete  entrainment  model.  Here  the  second  term  in 
eqn.  (28)  gives  the  growth  rate  of  Ah  and  thus  a  prognostic 
estimate  for  a. 

More  physics  may  be  added  by  revising  the  We  expression  to 
include  the  small  TKE  storage  rate  term,  Cslw,2/ih,  from  eqn.  (16). 
However,  our  present  focus  is  amending  w,,  rather  than  refining 
the  particulars  of  bulk  entrainment  models;  so  we  choose  to  rest 
with  the  above  diagnostic  and  prognostic  formulations  for  a. 

In  either  case,  hysteresis  in  a  versus  We/w*  is  evident,  and  like 
Kamada  (1988b,  figs.  1  and  2),  depends  on  changes  in  a0  and  z;. 
Both  models  include  shear  as  well  as  heating  due  to  entrainment, 
and  avoids  ad  hoc  surface  temperature  distributions,  or  fixed 
values  of  a,.  Also  note  that  eqns.  (16)  and  (21)  lets  us  specify 
a,  which,  together  with  eqn.  (22),  determines  D  from  eqn.  (4). 


BAROCLINIC  SHEAR 

If  the  geostrophic  wind  varies  linearly  with  height  according  to 
the  thermal  wind  relations,  u.  z  =  -(g/fT0)T  v,  and  vg>z  =  (g/fT0)T  x, 
then  from  eqn.  (23)  the  baroclinic  stress  profile  must  be 
parabolic,  rather  than  linear  as  in  the  barotropic  case. 

Wyngaard  (1985)  gave  the  following  simple  estimate  for  baroclinic 
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shear  based  on  the  thermal  wind, 


ub2  =  b,ug  2  +  b2vgz  *  vb,z  =  b;ug,z  +  bivg,z  /  (29a,  b) 

where  b,  =  1/ (m2  +  1),  b2  =  m/(m2  +  1),  m  =  fzt/w,,  and  f  is  the 
Coriolis  forcing.  If  so,  the  stress  can  be  written  as 


-u. '  w ' 


u,2  + 


b2z  (gz  p  f  T  y  +  2  fT  p  ) 
2  p  fT 


So  the  baroclinic  shear  contribution  becomes 


(30) 


h  _  gz,(b2Tx  -  b,Tv)  (b;gpz,:Tv  +  6pTu,2) 

1  -u.'w'u^az  =  - : - : -  , 

0  6  pfT-  (31) 

where  terms  in  b,  -  1  =  -m2/ (m2  +  1)  were  considered  neglectably 
small.  The  baroclinic  integral  is  taken  over  the  whole  vertical 
range  0  to  z,  rather  than  -L  to  Z|  because  the  SL  shear  analysis 
only  accounts  for  barotropic  effects  based  on  SL  similarity. 
Scaling  analysis  of  eqn.  (31)  shows  that  the  squared  and  cubic 
terms  in  z,  are  small.  If  so,  to  first  order,  the  baroclinic 
shear  contribution  from  eqn.  (31)  reduces  to, 


H 

f 


0 


-u, '  w  '  u,  tbz 


(z,/zb)w,3 


(32) 


where  we  can  define  a  baroclinic  length  scale  as, 

fT  w,3  w.3 

zh  =  - - -  =  - - -  (33) 

u,-g(b:Tx  -  b,Tt>)  u,2  ubz 

Usually,  zh  >>  zr  But  in  coastal  areas  mesoscale  temperature  and 
pressure  gradients  are  often  large.  For  a  coastal  z,  of  300  m,  a 
3  °C/100  km  horizontal  cross-wind  temperature  gradient,  We  near 
zero,  f  =  1 0‘4  sec1,  w,  —  0.8  m  s'1,  and  u,  =  0.4  m  s'1,  then  zb  = 

330  m,  and  w,s  may  be  enhanced  baroclinically  by  ~  50%. 
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SUMMARY 


In  summary,  we  account  for  surface  layer  (SL) ,  entrainment  zone 
(EZ) ,  and  baroclinic  (BC)  shear.  Each  source  is  given  in  terms 
of  a  dimensionless  ratio  of  the  relevant  length  scales  z0,  L,  ih, 
z,,  and  zb.  When  summed  and  multiplied  by  w, ,  this  results  in  a 
turbulent  velocity  scale  which  includes  both  buoyancy  and  shear 
generation  of  TKE,  namely, 


w,s  m 


n  -  (1.25  -  Ln  R)L/z1  +  O.SAh/Zj  +  zjzb 

n 


1/3 

w. 


(34) 


Here  n 


1  -  a,  (1  +  a) 
2 


al 


We 


WeA0/w'©'o  , 

(Cc  -  cd) 

- - -  w,3/Zi 

6a9  -  AUf/Ah 


the  Z;/surface  heat  flux  ratio, 


the  ML  entrainment  rate, 


a©  =  w'0'o/We  -  ruz,(l  +  q) 

2 . 3  (au2  +a  v2) 


a  s  Ah/z,  = 


I  ru  r 

L 


the  EZ  0  jump, 
1/2 


—  ,  the  EZ/BL  depth  ratio, 

2 


au  =  We(u«/fz,)2  ,  av  =  u«:/fz, 

R  =  (1  -  150z0/L)  1/3  -  1 


the  EZ  velocity  jump, 


an  SL  shear  generation  term, 


the  baroclinic  length  scale. 
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